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Abstract 



' We consider non-baclctraclcing random walk (NBW) in the nearest-neighbor setting on the Z''-lattice 

^ . and on tori. We evaluate the eigensystcni of the m x m-dimensional transition matrix of NBW where 

^ ' m denote the degree of the graph. We use its eigensystem to show a functional central limit theorem 

. for NBW on Z'' and to obtain estimates on the convergence towards the stationary distribution for 

NBW on the torus. 

(N 

1 Introduction 

Ph . The non-backtracking walk (NBW) is a simple random walk that is conditioned not to jump back along 
r~| ' the edge it has just traversed. NBW can be viewed as a Markov chain on the set of directed edges, where 

■ the edge has the interpretation of being the last edge traversed by the NBW, but we will not rely on 
this interpretation. We study NBWs on Z'^ in the nearest-neighbor setting, and NBWs on tori in various 
settings, derive the transition matrix for the NBW and analyse its eigensystem in Fourier space. We use 
this to study asymptotic properties of the NBW, such as its Green's function and a functional central 

J> . limit theorem (CLT) on Z'^, and its convergence towards the stationary distribution on the torus. In 
I particular, our analysis allows us to give an explicit formula for the Fourier transform of the number of 
f-f^ ' n-step NBWs traversing fixed positions at given times. We use this to prove that the finite-dimensional 

■ distributions of the NBW displacements, after diffusive rescaling, converge to those of Brownian motion. 
^SJ I By an appropriate tightness argument, this proves a functional CLT. We further evaluate the Fourier 

transform of the NBW n-step transition probabilities on the torus to identify when the NBW transition 
probabilities are close to the stationary distribution. Our paper is inspired by the study of various 
high-dimensional statistical mechanical models. For example, our derivation allows us to give detailed 
estimates of the probability that NBW on a torus is at a given vertex, a fact used in the analysis of 
^ ■ hypercube percolation in [6J. 

d . NBWs have been investigated on finite graphs in [3], in particular expanders, where it was shown that 

NBWs mix faster than ordinary random walks. See also |4j, where a Poisson limit was proved for the 
number of visits of an n-step NBW to a vertex in an r-regular graph of size n. In the nearest-neighbor 
setting on Z*^, NBWs were investigated in [SJ Section 5.3], where an explicit expression of its Green's 
function and many related properties are derived (see also [131 Exercise 3.8]). Noonan [10] investigates 
the generating function of NBWs, and his results also apply to walks that avoid their last 7 previous 
positions (i.e., with memory up to 8), and were used to improve the known upper bounds on the SAW 
connective constant. See [12j for an extension up to memory 22 for d = 2, further improving the upper 
bound on the SAW connective constant. The methods in [1^ \n\ allow to compute the generating 
function of the number of memory-A; SAWs for appropriate values of k, but do not investigate the number 
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of memory-/c SAWs ending in a particular position in U^. Finally, [TT] studies relations between the 
exponential growth of the transition probabilities of NBW and non-amenability of the underlying graph. 
For interesting connections to zeta functions on graphs, which allow one to study the number of NBWs 
of arbitrary length ending in the starting point, we refer the reader to [7]. 

This paper is oganized as follows. In Section [2] we investigate NBW on and in Section [3] we study 
NBW on tori. 



2 Non-backtracking random walk on Z 
2.1 The setting 

An n-step nearest-neighbor simple random walk (SRW) on 'Z'^ is an ordered (n+l)-tuple u = {ujq, uji,uj2, . . . , ojn), 
with LOi ^ Z and \\LOi — = 1, where ||x||i = Yli=i We always take (xIq = (0, 0, . . . , 0). The step 

distribution of SRW is given by 

D{^) = ^l{||x-|U=i}, (2.1) 

where 1^ is the indicator of the event A. If an n-step SRW u additionally satisfies that Ui ^ uji+2, then 
we call the walk a non-backtracking walk (NBW). As the problem of NBW is trivial for d = 1, we always 
assume that d >2. For the NBW we also count walks conditioned not to take their first step in a certain 
direction l. We exclusively use the Greek letters l and k for values in {—d, —d-\-l, . . . , —1, 1,2, ... ,d} and 
denote by G Z'^ the unit vector in direction l, i.e. (ej^ = sign(i)(5|^| (beware of the minus sign when i 
is negative, which is somewhat different from the usual choice of a unit vector). Let 6„(x) be the number 
of n-step NBWs with a;„ = x, and b'^{x) the number of n-step NBWs uj with Un = x and ui ^ e^. For 
n > 1, the following relations between these objects hold: 

E K-ii^-e.), (2.2) 
Le{±i,...,±d} 

b'ni^) + b-^,{x-e,) v., (2.3) 

E K^ii^-e.). (2.4) 
Ke{±i,...,±d}\{L} 

We analyse 6„ and b'^^ using Fourier theory. For an absolutely summable function / : Z'^ i— )• C, we define 
its Fourier transform by 

/>) = E/(^>''" (A;G[-vr,^]'^), (2.5) 

where k ■ x = ksXs and i denotes the imaginary unit. We use k exclusively to denote values in the 

Fourier dual space [— vr, tt]'^. For f,g: Z"^ i— )• C we denote their convolution hy f -kg, i.e., 

{f*g){x)=Y,fiy)9{x-y), (2.6) 

y&Z 
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bn{x) = 
bn{x) = 



and note that the Fourier transform oi f -k g is given by fg. Applying the Fourier transformation to 
dUI-dia]) yields 

1 

D{k) = -Y^cosih), (2.7) 

4=1 

bnik) = Yl k-iik)e"'% (2.8) 
te{±i,...,±d} 

bn{k) = b'^{k) + b-^^{k)e'^^ Vi, (2.9) 

Kik) = E K-ii^y'''- (2-10) 

Ke{±i,...,±4\{t} 

In our further analysis, we use C^'^-valued and C^'' x C^'^-valued functions. For a clear distinction between 
scalar-, vector- and matrix-valued quantities, we always write C^'^-valued functions with a vector arrow 
(e.g. v) and matrix-valued functions with bold capital letters (e.g. M). We do not use {1, 2, ... , 2d} for 
the indices of the elements of a vector or a matrix, but use {—d, —d+ 1, . . . , —1, 1, 2, . . . , c?} instead. Here, 
for a negative index l G {—d, —d + 1, . . . , —1} and a vector k G [— vr, vr]'^, we define k^ := — 

We denote the identity matrix by I G £'idx2d ^^^^ ^^le all-one vector by 1 = (1, 1, . . . , 1)"^ G C^'^. 
Moreover we define the matrices C, J G C^'^^^'^ by (C)^^^ = 1 and (J)^,^ = To characterize the 

displacement of a step in a direction i, we define the diagonal matrix D(/c) with the entries (D(A;))t^t = e''^'-. 

We define the vector bn{k) with entries ( bn{k) I = b'^{k). Then, we can rewrite ()2.8p - ()2.10p as 



bn{k) = T^t){-k)bn-l{k), (2.11) 

bn{k)l = bn{k) + -b{k)3bn-i{k), (2.12) 
bn{k) = {C-3)-b{-k)K-i{k) = ({C-3)-b{-k)Yl (2.13) 



We define the transition matrix 

A(A;) = (C - J)D(-A:), (2.14) 

so that (A(A:))t^K = e~*'''(l — S^-^)- With this notation in hand, we are ready to identify the NBW 
Green's function. 

2.2 The Green's function 

We start by deriving a formula for the NBW Green's function using the relations in (j2.1ip - (j2.1ip . While 
these results are not novel, the analysis presented here is efficient and simple. We define the NBW Green's 
function as the generating function of and b'^: 

oo oo 

4(x) := J2^bn{k)z^, B^^x) ■.= Y,K{k)z^ (2.15) 

ra=0 n=0 

B,{kf := {Bl{k),B-\k),Bl{k),...,B-\k)), (2.16) 
where denotes the transpose of the vector y G W^. By (|2.1ip - (|2.12p . 

B^{k) = 1 + zl^t){-k)B^{k), (2.17) 

Bz{k)l = Bz{k) + z'b{k)3B^{k) B^{k) = \l + zt){k)3] ^ lBz{k), (2.18) 

Using D(A;)JD(A;)J = I, it is easy to check that 

I + zD(fc)J ^ = - zD(A;)j) , (2.19) 
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and we use ([2T7ll - ([2J9]) to obtain 



^-(^) = T T^= , f (2.20) 



l-zl^B{-k) I + zD(A;)J 

Note in particular that 



f l + {2d-l)z'^ -2dzD{k) 



1-^2 



^.(^) = ^^(^^^C,.(.), (2.21) 
where = 2dz/{l + {2d - l)z'^) and C^{k) = 1/[1 - ^-D(/c)] is the SRW Green's function. By <^M), 



Mk) = I-zn{k)J lB,{k) 



I-z^ik) 
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l-z'l l + {2d-l)z^ -2dzD{k) 

so that 

B'(k) = . . (2.22) 

1 + i2d - l)z^ - 2dzD{k) 

2.3 The transition matrix 

The eigensystem. We start by evaluating the transition matrix (j2.14p by characterizing its eigenvalues 
and eigenvectors: 

Lemma 2.1 (Dominant eigenvalues). For d> 2 and k G {—tt,'k)'^, let A{k) be the matrix given in (|2.14|) . 
Then 



A±i = A±i(A;) = dD{k) ± ^{dD{k)f - {2d - 1) (2.23) 

are eigenvalues of A.{k). For k ^ (0, 0, . . . , 0)"^, the right eigenvectors {7**^' to the eigenvalue X±i are 
given by 

i;(±i) = A±ir±D(fc)l. (2.24) 
For fc = (0, 0, . . . , 0)-^, the eigenvectors are given by ^T*^'(0) = (2d — 2)1 and v^~^^ := (1, —1, 0, 0, 0, . . . , 0) G 

As we will see below, Lemma 1 2. II vields the two most important eigenvalues. When A-|- = A_, which 
occurs when {dD{k))^ — {2d — 1) =0, it turns out that has geometric multiplicity 1, and that 1 
is a generalized eigenvector satisfying A(A;)1 = tT*^' + A+l. We continue by computing the remaining 
eigenvalues and -vectors: 

Lemma 2.2 (Simple eigenvalues). For d > 2, k £ (— 7r,7r)'^, let A{k) be the matrix given in ()2.14p and 
el € C'^'^ the ith unit vector, i.e., (ej^ = J^.k for k G {±1, . . . , ±d}. 
For i G {2,3, . . .,d}, let 

^0 = (i + e'*^')(e''''ei + e_i)-(l + e''=i)(e^^^e; + e_J VA: G [-vr, ^]'^ 

= (l-e^^O(e''''ei-e_i)-(l-e^'=i)(e^^'e, -e_J ^k e [-Tr^Trf if e'''' ^ 1, 

Then, iJ'*'' is an eigenvector of A{k) to the eigenvalue =Fl- Both eigenvalues have a geometrical multi- 
plicity of d — 1 for all k. 

Lemmas I2.1H2.2I identifv a collection of 2d independent eigenvectors, and thus the complete eigensys- 
tem, of A{k). Now we prove these two lemmas: 
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Proof of Lemma [2A[ Let A G {Ai, A_i} and -u = Al — D(A;)1. The values Ai and A_i are the solutions of 
the quadratic equation 

P = 2d\b{k) - {2d - 1). (2.25) 
Using CD(-A;)1 = 2dD{k)l and JD(-A;) = D(A;)J, we compute 

A{k)v = {C-3)t){-k){XI-'b{k))l= (2d\D{k)I-X'D{k)-{2d-l)I^l 
'Pl- AD(A:)) l = Xv. 



This proves that iT is a eigenvector of A{k) corresponding to the eigenvalue A for all A; 7^ and also for 
the case of z. = 1 for k = 0. For A; = (0, . . . , 0) we note that A_i(0) = 1 and see that 

A{oy-'\0) = (C - 3y-'\0) = -3v^-'\0) = v^-'\0). 

□ 

Proof of Lemma{27^ For i G {1, 2, . . . , d}, the vectors 

= e'^^'^el + e_t and n^"'^ = e^'^'e^-C-, 
are eigenvectors of JD(— A;), where -5*='=''' is associated to the eigenvalue ±1. For l G {2, 3, . . . , d}, we define 
_ ^(1) _ ^(0 ^-0 = ^{-D Y n^T'' - n^"" ^ 4"^). (2.26) 



By construction, i;*'' and v'- are also eigenvalues of JD(— A;). For CD(— A:), we compute that 

CD(-A;)n(^' = C{e, + e'^^'e.,) = ^n^l, CD(-A;)n<-" = C{e, - e'^^'e.,) = - ^nlr"!, 



so that 



Knowing this, it follows that 



CD(-A;)iTW = CD(-A:)f;<-'' = 0. 

A{ky''> = (C - 3)'b{-k)v^'^ = -v^'\ (2.27) 
A(A:)i;<-'' = (C - J)D(-A;)iTW = ^-o_ (2.28) 

By ()2.27p . v^'^ is an eigenvector for all k. Since the set of vectors (?7''^)t=2,3,...,d is linearly independent we 
know that the eigenvalue —1 has geometric multiplicity d — 1. 

From (|2.28p . we conclude the existence of d — 1 linear independent eigenvalue for 1 only when e''^' 7^ 1 
for all i G {1, 2, . . . , d}. To prove that the eigenvalue 1 has geometric multiplicity d — 1 for all k, we show 
how to choose d — 1 linear independent eigenvectors when e^'^" = 1 for a. k ^ {1,2, . . . ,d}. For this, let Si 
be the set of all k G {1, 2, . . . ,d} with A;^ = and S2 the set of all k G {1, 2, . . . ,d} with A:^ 7^ 0. Let si 
and S2 be the number of elements in 5i and 82- Then S1 + S2 = d. For all t G ^i, we define v'-''^ = e^ — e_t. 

If si = d, then A; = 0. In Lemma l2.H we define for this case A2 = 1 and = ei — e_i. Then 
|^(-i)^^-2)^ . . . is a set of independent eigenvectors of A(A;) to the eigenvalue 1. 

If si < d, then let p be the smallest number in 5*2 and define 

for all i G 5*2 \ {t}. Then it is easy to verify that the vectors (i^'~'')t=2,3,...,d ^-re linearly independent and 
are eigenvectors of A (A;) with eigenvalue 1. □ 
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We use the eigensystem of the matrix A(A;) to identify bn{k) and bn{k): 

Lemma 2.3 (NBW characterization). Let d > 2, n > 1 and k G (— 7r,7r'^) such that Xi{k) ^ A_i(fc). 
Then, 

. , D(A;)(A^^(A;) - A'^i(fc)) + A"t1(A;) - Ari(fc) , , 

bJk) = 2d ^ 1^ ^ ."^^ ^' , — 2.29 

Ai(fc)-A_i(A:) ^ ^ 

bJk) = .r^i^D fc 1 - \ ^ ^ . ^ ^ 1. 2.30 

Ai(fe)-A_i(A;) Ai(fe)-A_i(A:) 

W^/ien Ai(A;) = A_i(A;), 

bn{k) = 2d[{n - l)Ai(A;)"-2 + D{k)nXi{k)''-^], hn{k) = [{n + l)Ai(A:)" + nAi(A;)"-iD(fc)]l. (2.31) 
Clealry, one can reprove (|2.20p and (|2.22p using Lemma l2.3i 
Proof. For Ai(A;) ^ A_i(/c), we define 

a{k) ^ ^ 



\i{k)-\.i{k) 2^J{dD{k)y - {2d -I) 
We can write 

1 = a{k)(^Xi{k)I-\^i{k)I + B{k)-i){k)'^l = a{k)v^'\k)-a{k)v^-''>{k). 

Using ()2.13p and the fact that v'-'^^^{k) are eigenvectors of A.{k) with eigenvalue A-|-(A;), we obtain 

bn{k) = A{kri = a{k)(^X'^{k)if^''>{k)-X1^{k)v^-''>{k)'^ , (2.32) 

which proves ()2.30p . Combining ()2.1ip and ()2.32p gives 

h^{k) = a{k)l^t){-k) (^XUky\k) - X1^{ky-'\k)y (2.33) 

Inserting the definition of v'-''^^^{k) gives (j2.29p . The proof of ()2.3ip is similar, now using that A(A;)1 = 
u'^' + A+1, which implies that A(/c)"l = nA+(/fc)"-^{;'i) + A+(/c)'"l. □ 

2.4 Central limit theorem 

This section is devoted to the proof of a functional central limit theorem for the NBW. The uniform 
measures on n-step NBWs form a consistent family of measures, so that there is a unique law that 
describes them as a process. We let to = {u}r)r>o be distributed according to this law. For a NBW oj and 
t > 0, we define 

Xn{t) = (2.34) 



where [xj denotes the integer part of x G M. 

Theorem 2.4 (Functional central limit theorem). The processes (Xn{t))t>o converge weakly to {B{t))t>o, 
where {B{t))t>o is a Brownian motion with covariance matrix!/ [d — 1). 

Our proof of Theorem 12.41 is organized as follows. We use Lemma [2. 3 1 to prove a CLT for the endpoint 
in Lemma 12.51 We then prove the convergence of the finite-dimensional distributions to the Gaussian 
distribution (see Lemma 12. 6p , followed by a proof of tightness (see Lemma 12. 8p . This implies Theorem 
12.41 see e.g. [51 Theorem 15.6]. We now fill in the details. 
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Lemma 2.5 (Central limit theorem). Let d>2 and k G [— vr, vr]'^. Then, 

lim Kle'^-^"/^] = e-myi2<i-2)^ .2.35) 

n—^oo 

where \\k\\2 = Yli=i Euclidean norm of k. 

Lemma 12.51 implies that the distribution of the endpoint of an n-step NBW converges in distribution 
to a normal distribution with mean zero and covariance matrix [d — 1)~^I. 

Proof. We can rewrite the expectation as 

]gr ifc-aj„/v^i ^ ^n{x) ik.^/^ ^ bn{k/^/n) 

L ^ ^6„(0) 2di2d-ir- ^ ■ ' 

As n — )• oo, we can assume without loss of generality that k is small and therefore that Ai(A;) > A_i(A;). 
Then we can use (|2.29p to compute the limit of ()2.36p . For i = 1,-1, we compute 

(2d-i)»-i "^^/^^ Yd • ('-^^^ 

We first consider the case l = \. The coefficient a{k) as well as v'^^^{k) are continuous in a neighborhood 
of 0, so we can directly compute 

rU{-k/VE)v^^k/VE) 1^D(0)^-W(0) 

Yd = Yd = 

Further, Ai(A;) is differentiable in k, so we can Taylor expand Ai(A;) at to obtain 

x,[k) = 2d-i-^^\\k\\i + oml)- 



Therefore, 



n^oo (2d - l)"-i n->oo V2d-2n \ J J 

so that the limit (piBT]) for t = 1 is given by e^ad^"''"^. 

We next consider the case i = —1, for which we use that k i— >• A_i(/c) is continuous and A_i(0) = 1, 
Therefore, for d > 2, 

lim -\ ^ , ; = 0. 2.39 

n^oo {2d - !)"-! ^ ^ 

The second factor in (|2.37p can easily be bounded uniformly for small k, so that for t = — 1 the limit of 
(p:3711 is zero. □ 

Lemma 2.6 (Convergence of finite-dimensional distributions). For d > 2 and N > 0, let = to < ti < 
t2 < ■ ■ ■ <tN = ^ and fcM G (-vr, vr]'^ for r = 1, . . . , N . Then, 

lijj, E[gi(E^=lfc('->('^LtrnJ-'^LV-inj)/v^)] = q- Jlr^l W^^''^ WUtr-tr -l) / {2d-2) _ (2.40) 



n— >oo 
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Proof. As we take the limit n — ^ oo, without loss of generality we can assume that r]r{n) := [trn\ — 
[tr-in\ > 1, tr > tr-i and each kn^ := k'''''>/^/n is so small that Xi{k) > A_i(/c) for r = 1,...,A^ and 
n G N. Let Wn be the set of all n-step NBW. For any function /: Z''^^ H> C, we know that 



IE[/(w[tl„J, • • • ,W[f^nj)] = X] fi^ltin\, ■ ■ ■ ,l^tNn\) 



TTTs Yl Yl f{xi,...,XN) Yl 



bn{0) 



i=l N 



(2.41) 



Let hn^{x) be the number of n-step NBW oo with wi 7^ e^, a;„_i = x + e^ and = x. We define the matrix 
B„(A;) with entries (B„(A;))t,K = bn'^{k). By a relation similar to ()2.10p . we conclude that B„(A;) = A(A;)". 
We fix N points xi, . . . ,xj\f £ Z, then the number of NBW lu with uj^^^n] = fo^ alH = 1, 2, . . . , is 
given by 



in {±l,...,±d} 



r=2,...,Af 



(2.42) 



We insert /(xi, . . . ,xn) = e'^'-'^" (xr-xr-i) g^^^ obtain 



1 



^n(O) 

in {±l,...,±d} 



N 



r=2 



^rD(-A;W)A(fcW)^i(")-i nB,^(„)(A;W)r 

On(0) 



r=2 
TV 



^FD(-A:«)A(A:«)^iW-inA(A:r)''^(")r. 

On(0) 



r=2 



(2.43) 



To proceed, we define, for r G {1, 2, . . . , A^}, 



Kir) 



FD(-A;(,^))A(fc(,i')''i(")-^ • • • A(A;(,^) )''-(") 



K{0) = l™(-fc«). 



Z]r=i ^»'('T') < ^T^^, by construction of hnir), for all v G C^"^ 



|/in(T)i;| < E 



(2.44) 
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where we write ||i'||oo '■= max^ \vi,\. Therefore, we can rewrite 



2d -I 2d -2 '^-^ ' y 2d -2 

U^)^{^]\{C^Y'^' (2.45) 



„^ „ ' '2T- 

r=l i=r+l 



To compute the Umit of (j2.43p . we show that (|2.46p and (j2.47p converge to zero and identify the limit 
of ([235]). From ([Tli]) . it follows that we can bound ([2^ by I - v'^^\k''n^) / {2d - 2) and we know that 
{/'^^(fcn^') {2d — 2)1. Thereby we conclude that (|2.47p converges to 0. 

To compute the limit of (IMD we note that Ai(A;) <2d-l and that v'^'+^\k'il+^^) - v^^'^k^n'') 0, so 
that 

hm IdZaSDI < hm ^fJ- (^-<^)(A;«) - ^^\k^^)) 

r=l 



< lim — -^|b-(^)(fc«)-T;«(fc<^))IU = o. 

rn-oo Zd — 2 ^ — ' 
r=l 

To compute the limit of (|2.45p . we use that 



r .^^^*'HfcnO 2d-l i^D(-A;)(Ai(A;^^')r+ D(A;^^')r) 2d - 1 ^ 

lim n„{U) ; 7— = hm ; — ; ^ -— = 1, 

^ 2d -2 X^(ki!^) n^cx. 2d(2d-2) Ai(e') 



and, for each r G {1, . . . , A^}, 



r]r{n) 



lim =6 

n->oo \ 2a — 1 / 

Combining this yields that the limit of ()2.45p is the right-hand side of ()2.40p . which completes the proof 
of Lemma 12.51 □ 



To prove tightness we make use of the following lemma, which computes the second moment of the 
end-point of NBW. The leading order in this result was alternatively proved in [9j (5.3.11)] using residues. 

Lemma 2.7 (The second moment). For d > 2 and n G N, 

2, d M-l d 



lE[||^^n||2] = 3 7^ + + 



d-1 2{d-lY 2(d- l)2(2(i- l)"-2' 
Proof. We define the differential operator := f M") ^"^^ t^aX 



dki 

nWn\\l] = J—Y.hr.{x)\\x\\l = -^V26„(0). 

6n(0) 6„(0) 
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To compute the second derivative in a neighborhood of the origin, we recaU ()2.29p . By Lemma 12.31 ^oi k 
such that X+{k) / X~{k), 



(2.48) 



'^£{-1.1} y (.dD{k))'^ - {2d - I) ae{-i,i} 



Since, A+(A;) ^ ^-ik) for k smah, (j2.48|) holds in particular for k small. A straightforward computation 
gives 



k=0 



da V^D(O) 



(dL'(0))2 - (2d - 1) 



d^D{0) 



-da 



(dZ)(0))2 - {2d-l) 
d{{l + a)d - 1) nAJJ-i(O) - (n - l)A^-2(0) 



A-(0)i^(0)-Ari(0) 



{db{fd)f - {2d -I) 



da 



d-1 



(d-l) 



d-1 

2 Aa(0)-l ' 

(d-l)2 



so that 



We arrive at 



vV(fc) 

VVi(A^) 



fc=0 



d-l d-1 



fc=o d — 1 



(d-1) 



[2d(d- l)n + (4d- 1)] 
, 2d - 1 



V'6„(0) 



2d(2d 



_ l)n-l 



-vV(fc) 



d 4d- 1 

-?^+ TT^ TTTT + 



k=0 



d-1 2(d-l)2 2(d- l)2(2d- l)'^-2' 



□ 



Lemma 2.8 (Tightness). For d > 2 and < ti < t2 < < 1, there exists a K > 0, such that, for all 
n>l, 

E[\\Xn{t2) - Xn{tl)\\l\\Xn{t3) - Xn{t2)\\l] < K{t2 - ti){t^ - ta). 

Proof. We use the same notation as in the proof of Lemma 12. 6i In ()2.42p , we have seen how to describe 
the number of NBWs that visit a number of fixed points. This time we forget about the non-backtracking 
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constraint between two subsequent NBWs to upper bound 



- TTT^ X] \\X2- Xi\\l\\x3- X2\\lb^^(n){xi)Y[br,^(n)iXr - Xr-l) 

— X] Ik2||2l|2;3||2n 



xi,X2,X3,X4&Z'^ r—1 



^m(n)^V4,(n) , , nO r / \ii ii2 
1 2^ \^(n)\X2)\\x2\\2 2^ 0,,3(n)(^3)||a;3||2 

IE[||tj^2{n)ll2]E[l|wr,3(n) nil- 
Applying Lemma 12.71 completes the proof. □ 

2.5 Extension to non-nearest-neighbor setting 

In this section, we extend the analysis of NBW on Z"^ to other bond sets. We start by introducing the 
bond sets that we consider. We let IB C Z'^ x Z'^ be a translation invariant collection of bonds. Let 
Vo = {x: {0, x} G B} denote the set of endpoints of bonds containing the origin and write m = |Vo|- We 
assume that Vq, and that Vq is symmetric, i.e., —x S Vq for every x G Vq. Thus, m is even. We 
define the simple random walk step distribution by 

D{x) = ^l{.eVo}- (2.49) 

Define the matrices C, J S C*"^™ by (C)^^^ = 1 and {J)x,y = ^x-y^ and let the diagonal matrix D(fe) 
have entries {D{k))x,x = e''^'^, where x,y G Bq- Then, we define the matrix A.{k) of size m x m by 

A(A;) = (C - J)D(-A;). (2.50) 

With this definition at hand, we see that (|2.11|) - (|2.13|) remain to hold. As a result, also Lemmas 12. lH2.2t 
whose proof only depends on ()2.1ip - ()2.13p . continue to hold when we replace each occurrence of 2d by m. 
Since the proof of Lemma 12.31 in turn, only depends on Lemmas I2.1H2.2| also it extends to this setting, 
so that, for example 

^±{k) = F±(D{k);'m), where F±{x; m) = — (jnx it -s/ {mxY — 4(m — 1)^ , (2-51) 

and, when Ai(A;) / A_i(A;), 

. . . ^ m mmk) - A!!i(fc)) + fe^(fc) - XrHk)) 

Naturally, Theorem 12 . 4 1 needs to be adapted, and now reads that the processes {Xn{t))t>o converge weakly 
to a Brownian motion with covariance matrix M of size d x d, where, for t, n £ {if; • • • > ^d}, we define 



M, 



d^Xiik) 



XiiOr\ (2.53) 

fe=0 



dk.dkf. 

We next compute M explicitly in terms of the covariance matrix of the transition kernel D. We compute 
that F+(l;m) = m/2 + (m - 2)/2 = m - 1, and 



m?x 



F|(rc;m) = m/2 H , so that F+(l; m) = m(m - l)/(m - 2). (2.54) 

2-y/ {mxY — 4(m — 1) 
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By symmetry, the odd derivatives of D{k) are zero, so that a Taylor expansion yields 

i)(A;) = 1 - ^A;^HA; + 0(||A:||^), (2.55) 

where, for i, k £ {1, . . . , d}, 

li,,n = Y,x,x,D{x) (2.56) 

X 

denotes the covariance matrix of SRW. As a result, 

F' (1- m) 

M = U-^ — f = Hm/(m-2). (2.57) 
F+(l;m) > ^ ' 

In the nearest-neighbor case, m = 2d and H = 1/d, so that we retrieve the result in Theorem 12.41 

3 NBW on tori 

In this section, we extend the results in Section [2] to NBWs on tori. In Section [3 . II and [3 . 2 1 we investigate 
NBW on a torus of width r > 2, and in Section 13.31 we investigate NBW on the hypercube, for which 
r = 2. The study of random walks on various finite transitive graphs has attracted considerable attention. 
See e.g., [H] for a recent book on the subject, and [2j for a book in preparation. Here we restrict ourselves 
to NBWs on tori. 

3.1 Setting 

For d>2 and r > 3, we denote by T = T.^^^ = (Z/rZ)'^ the discrete d-dimensional torus with side length 
r. The torus has periodic boundaries, i.e., we identify two points x,y £ Tj. ^ if Xi mod r = yi mod r for 
all i = 1, . . . , d where mod denote the modulus. We define the Fourier dual torus of T as 

r I 

so that each component of G T* ^ is between — vr and vr. The Fourier transform of / : — t- C is defined 

by 

f{k) = J2 /(^)«*'' ^ ^ ^M- 

As in Section [21 we define an n-step random walk on T^^^ to be an ordered tuple uj = (wq, . . . ,w„), with 
G 1'r-,d and LJi — Wj+i G Vo, where we recall that Vq = {x: {0, x} G B}, and B is the translationally 
invariant bond set on which our random walks moves. We always assume that Vqi and that Vq is 
symmetric, i.e., if x G Vq, then also —x G Vq. Further, we always assume that uq = (0, . . . ,0). The 
simple random walk step distribution is given by 

D{x) = -%.ev4 and D{k) = - Yl (3.2) 

If an n-step random walk on T^^d additionally satisfies 7^ then we call the walk a non-backtracking 
walk (NBW) on T^. ^^. Let 6ri(x) be the number of n-step NBWs with a;„ = x. Further, let h':J^x^ be the 
number of n-step NBWs a; with a;„ = x and wi 7^ . 

In this setting, WG ca,n express h^S^x) for NBW on T^^. in terms of NBW on Indeed, identify Ty.^^ 
with {0, . . . , r - l}'^ C Z'^, and also identify Vq = {x : {0, x} G B} as a subset of {0, . . . , r - 1}'^ C Z^. 



(3.1) 
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Define Vg = Vq U (~Vo) (which, by construction, are disjoint subsets of and define the random walk 
step distribution D-^ix) by the uniform distribution on Vq. Then, for x G T^^d, 

K{x) = hl{y), (3.3) 

y: y^x 

where, for y G Z'^ and x G T^^d, we say that x ~ y when x = y mod r, and h^{y) denotes the number of 
n-step NBWs on with step distribution Dj,{x). As a result, 6„(A;) for NBW on T^^d is equal to h\{k) 
for every /c G T*^ as defined in (|3.ip . Therefore, we can use most results for NBW on Z"' to study NBW 
on Ty^d- We define the probability mass function of the endpoint of an n-step NBW by 

bn{x) bn{x) 

Pn{x) = ^ , . . = —f tt;^. (3.4 

L.y bn[y) rn[m - 1)" i 

In order to study the asymptotic behavior of NBW, we investigate \±{k) for k ^ 0. Our main result 
in this section is the following theorem: 

Theorem 3.1 (Pointwise bound on bn(k)). Let Yq be symmetric and satisfy Vq. Then, for n G N, 
NBW with steps in Vq satisfies 

\Pn{k)\ < {l/V^^^V\D{k)\y-\ (3.5) 

To prove Theorem 13. 1|. we start by investigating A-|-(A;) for k ^ 0. For this, we use Lemma [2. II to note 
that ^ 

^±ik) = F±{D{k);m), where F±{x;m) = — (^mx i: \l (mx)"^ — 4(m — 1)^ , (3.6) 

and where m denotes the degree of our graph. We bound X±{k) in the following lemma: 
Lemma 3.2 (Bounds on A-|-(A;)). For any k ^T*^, 

= \/m — 1 when {mD{k))'^ — 4(m — 1) < 0; 

\X+{k)\ < < (m - 1)[1 - (1 - b{k))m/{m - 2)] when D{k) > 0, {mb{k)f - 4(m - 1) > 0; (3.7) 
^ < 1 when b{k) < 0, {mb{k)f - A{m - 1) > 0. 

Proof. The function x i— )• F±{x; m) is real when (mx)^ — 4(m — 1) > 0, and complex when (mx)^ — 4(m — 
1) < 0. When {mxf' - 4(m - 1) < 0, 

|F±(x;m)p = m- 1, (3.8) 

so that \F±{x;m)\ = \Jm — 1. 

When (mx)^ — 4(m — 1) > 0, by the symmetry i<+(x;m) = F^[—x;m), we only need to investigate 
X G [0,1]. We start with F-{x;m), which clearly satisfies F-{x;m) > 0. Further, we can compute that 
F_(l; m) = 1, and 



m'^x m 



mx 



F'_{x-m) =m/2 = — 1 -^^^=^^^= < 0, (3.9) 

2y'(mx)^ — 4(m — 1) 2 L y'(mx)^ — 4(m — 1) ^ 

so that F^{x] m) < 1 for all x G [0, 1] for which (mx)^ — 4(m — 1) > 0. 
To bound F+(x;m), we use ()2.54p as well as 

2 4 2 

TP"/ ^ m X 

F^{x; m) = ^^^^2 _ 4(^ _ i) " 2((mx)2 - 4(m - l))3/2 ^^■^^> 
rr? 

2((mx)^ — 4(m — Ij)'^/^ 

2vr?{m — 1) 

= < 

((mx)2 - 4(m - l))3/2 
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As a result, a Taylor expansion yields 

F+(x;m) < F+(l;m) + (x - l)F\_{l;m) = (m - 1) + (x - l)m(m- l)/(m - 2) (3.11) 
= (m — 1)[1 — (1 — x)m/{m — 2)]. 

□ 

Proo/ of TheoremlMl By ([2111) and (f233D . 

|6n(A:)| < ||D(-A;)r||2||6„_i(A;)||2 < UllbH A(A;)||Sp^||r||2, (3.12) 

where we write ||M||op = sup ||Mx||2/||x||2 for the operator norm of the matrix M. We next use that 
A(A;) has eigenvalues A+(/c),A_(fc) and ±1 by Lemmas [2?T][2i2l so that 

||A(A:)||oP = |A+(A;)| V|A_(A;)| VI, (3.13) 

where we use that for finite-dimensional matrices, the operator norm is equal to the maximal eigenvalue, 
and for x, y G M, we write (x V y) = max{x, y}. Thus, we arrive at 

\bn{k)\ < m(|A+(A:)| V |A_(A:)| V l)""\ (3.14) 



By Lemma |3.2|, and since m > 2 so that \/m — 1 > 1, 

MM < _ i)-i/2 V (1 - [1 - ^(A:)]^^^) < (m - l)-i/2 V |^(A:)|. (3.15) 
m — 1 ^ m — 2' 

Substitution into ()3.12p yields the claim. □ 
3.2 Asymptotics for NBW on the torus 

In this section, we study the convergence towards equilibrium of NBW on tori of width r > 3. We focus 
on two different examples. The first is random walk on products of complete graphs, where 

Vo = {x: 3!i e {1,... ,4 such that x-i / 0}. (3.16) 

Our second example is NBW on the nearest-neighbor torus. The reason why we study these cases 
separately is that NBW on products of complete graphs is aperiodic, while nearest-neighbor NBW is 
periodic. Therefore, the stationary distribution for NBW on products of complete graphs equals the 
uniform distribution on the torus, while for nearest-neighbor NBW, the parity of the position after n 
steps always equals that of n. In Section 13.31 we further study random walk on the m-dimensional 
hypercube. 

For any small ^ > 0, we write Tinix(^) for the ^-uniform mixing time of NBW, that is, 

rmix(4) = min<^ n: max < (1 + 4)^ h (3-17) 

where V = r'^ Is the volume of the torus. We start to investigate NBW on products of complete graphs: 
NBW on products of complete graphs. Our main result is as follows: 

Lemma 3.3. For every d > 1, r > 3, n > log ((r - 1)/[(1 + 0^'^ _ i])^ nbW on products of 

complete graphs satisfies that 

max \pn{x) - r-'^\ < (m - l)-(«-i)/2 + ^^-d. (3.18) 



In particular, for every e > 0, there exists Vq such that when r'^ >Vo, 

TrnUO < (1 + £)^^^M) log ((r- - 1)/[(1 + ^i/rf _ 1]) . 
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When ^ is quite small, we obtain that T,^ix{0 ^ (1 + log {d{r — 1)/^). 

Proof. The inverse Fourier transform on T^^d is given by 



r,d 



SO that 



Therefore, 



^"(^) = 4 E ^-(^>''"=^"'^"(0) + ^ E (3.19) 



To bound we rely on Theorem 13.11 and start by computing 



d r-l 



^W = ^^EE«^'"- (3-21) 
Since ki e ^ {- [^\ \^]}, we have that ^^^"^ e'^^^^' = for ki / 0. Therefore, 

j=i j=i j=i 

where o(A:) = X^f^x Ijfci^o} denotes the number of non-zero coordinates of k. By this observation, 
^{k{j)) = 1 — ^^J^^j^^ for any k[j) G T*^ for which a{k{j)) = j. Then, by Theorem 13.11 and the fact 
that there are {'f)ir — ly values of A; G T*^ for which a{k{j)) = j, 



d 

d\ „-d 



\pn{x) -r \ = r 



V \Dikimr-' (f) {r - ly (3.23) 

7 = 1 ^-^^ 



J 

d 



< _ l)-(n-l)/2 + ^-rf ^ J^^^ _ l)ie-n(n-l)/[d{r-l)] ^ 



where we use that |1 — ^i^^zrjyl < e~^-^^"~^^/['^(^~^)l for any j = 1, . . . ,d. Thus, 

\pn{x) - r-'^\ <{m- l)-("-i)/2 + r-'^[(l + (r - i)e-'''^n-i)/ld(r-i)]^d _ 

< (m- l)-("-i)/2 ^^r"*^, (3.24) 

when n > "^^ log ((r — 1)/[(1 + S,y^'^ — 1]). The result on T^i^iO follows immediately. □ 
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NBW on the nearest-neighbor torus. Our main result is as follows: 

Lemma 3.4. For every d > 1, r > 3 and n > log(2/[(l +^/2)^/'^ - 1])/[1 - cos(27r/r)], NBW on the 
d- dimensional nearest-neighbor torus satisfies that 



max \pn{x) - [1 + (-l)ll^lli+"]r-'* < {2d - 1)""/^ + ^r-". 

In particular, for every e > 0, there exists Vq such that whenever r'^ > Vq, 

TmUO <(! + £) log (2/[(l + C/2)i/^ - !])/[! - cos(2^/r)]. 
When r'^ is large and small, the above implies that T^ixiO < (1 + 2e)(r^/(27r^))(log {2d/S)). 
Proof. We adapt the proof of Lemma 13.31 to this setting. We have 



(3.25) 



(3.26) 



Further, D{-n) = -1, so that, by (ICTD bni^) = 6„(0)(-l)". Therefore, 



,(x)-[l + (-l)INIi+"]r-'^| = 1 



^n(^)gi/c-z 



< 



- y 



\K{k)\ 



<(^_l)(n-l)/2+ 1 ^ |^(fc)r^ 



In the nearest-neighbor case, D(k — tt) = —D(k), so we may restrict to k for which D{k) > 0. Let 

T:,d,+ = {A:GT;_,:^(fc)>0} 
denote the set of fc's for which -D(/c) > 0. Then, 



We use that 
so that 



|p„(x)-[l + (-l)ll^lli+"]r-'^| < (m-l)("-i)/2^— Yl D{kT~^. 
b{k) = !-[!- i){k)] < e-[i-^W], 



,(x)-[l + (-l)ll^'lli+"]r-^| <(m-l)("~i)/2^— Y 



feGT* : fc^O 



(m-l)('^-)/^ + -|[ 5: e-("- 



{n-l)[l-D{fc)] 



l)[l-i5(fc)] 



(3.27) 

(3.28) 

(3.29) 
(3.30) 



= l)(«-l)/2 + ^^(^ ^ g-(n-l)[l-cos(fc)]/d^'^_ 

We use that the dominant contributions to the sum over /c G T* ^ comes from k = Q and k = ib27r/r, so 
that 

g-("-i)[i-cos(fc)]/d _ ]^ _|_ 2e~("^^)[^^™^(2'^/'')l/'^)(l + o(l)), (3.31) 
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so that 



•n{x) - [1 + (-l)ll^lli+"]r-'^| < (m - l)("-i)/2 + _ + 2e-("-^)[^-™'^(2^/'')l/'^)(l + o(l))) - l' 

< (2(i-l)-("-^)/2^^r-'^ (3.32) 

when n > (log (2/[(l + C/2)^/'^ - 1])/[1 - cos(27r/r)]. □ 
3.3 NBW on the hypercube 

In this section we speciahze the results of Sections [2] and Sections 13. 1H3.2I to the hypercube T2^m = {0, 1}™. 
The results in this section are an important ingredient to the analysis of percolation on the hypercube in 
[6]. We start with some notation. It will be convenient to let the Fourier dual space of {0, 1}™" = {0, 1}"* 
be = {0, 1}™, so that the Fourier transform of a summable function / : {0, 1}™ — t- C is given by 

/>)= E /(^K"'"= E {k€{o,ir). 

xe{o,i}"^ xe{o,i}™ 

For k £ {0, 1}™', let a{k) be its number of non-zero entries. Then, the SRW step distribution on {0, 1}™ 
satisfies 

^ m 

b(k) = —Yi-^f' ='^-Mk)/m. (3.33) 
m ^-^ 

3 = 1 

The main result of this section is as follows: 
Theorem 3.5 (NBW on hypercube). For NBW on the hypercube {0, l}*", 

Pn{k) < (\b{k)\ V 1/Vm - l)" \ (3.34) 
Consequently, for every e > 0, there exists mo such that for all m > mo, 

TmUO < "^^^^'^ log {2m/ 0. (3.35) 

Random walks on hypercubes have attracted considerable attention, In fact, the bound on the uniform 
mixing time for NBW in discrete time closely matches the one for SRW in continuous time (see [U Lemma 
2.5(a)]). 

Proof. The bound in (j3.34p follows directly from Theorem l3.1l . We continue to investigate the convergence 
of the NBW transition probabilities x i— )• Pnix) to its quasi-stationary distribution, which is 2 x 2"™ when 
n and x have the same parity and otherwise. Here we say that n and x have the same parity when 
there exists an n-step path from to x. 

Lemma 3.6 (Convergence to equilibrium for NBW on hypercube). For NBW on {0, 1}™ and every 
n > d(log(i-hlogO/2, 

max \pn{x) - 2-'"[l + (-l)ll^lli+"]| < (m - l)-("-i)/2 + ^2""^. (3.36) 
xe{o,i}™ 

Consequently, for every e > 0, there exists mg such that for every m > mo, 

TmUC) < -^^^^^log([l +C/2)]V- - 1). 
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Proof. The inverse Fourier transform on the m-dimensional hypercube is given by 

fix) = 2-^ 

A:e{0,l}™ 

so that, using bn(l) = (— 1)"6„(0), 

6„(a;) = 2— ^n(A:)(-l)'-" = 2-'"[l + (-l)Wi+"]6„(0)+ ^ (3.37) 

fce{o,i}™ fce{o,i}™: fc^^i.o 

Substituting the bound ()3.34p in p.37p leads to 

m/2 



Pnix) - 2— [1 + I < 2 • 2— ( ) [(1 - 2j/mr~' + (m - 1)""^ 

i=i ^ ^ 

m/2 , , 

< _ i)-/2 + 2 • 2-"^ ^ ( j e-2^("-i)/™ 

j=i v.? / 



(3.38) 



< (m - l)-"/2 + 2 • 2-™[(l + e-^C"-!)/™)"^ _ l] 

< (m-l)-"/2+^2-™, 

when n > -flog ([1 + e/2]^/'" - 1). □ 
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